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ABSTRACT 

We model the sub-millimetre polarization patterns that are expected for filamentary clouds 
that are threaded by helical magnetic fields. We previously developed a three parameter model 
of such clouds (Fiege & Pudritz 2000a) , which are described by a concentration parameter C, 
and two flux to mass ratios T z and to specify the mass loading of the poloidal and toroidal 
field lines respectively. Our models provide a simple and purely geometric explanation for 
the well-known "polarization hole" effect, in which the sub-millimetre polarization percentage 
decreases toward the regions of peak intensity. This occurs because of a cancellation between 
contributions to the polarization from the "backbone" of poloidal flux along the filament's 
axis and its surrounding envelope, which is dominated by the toroidal field component. A 
systematic exploration of our parameter space allows us to classify the polarization patterns 
due to filaments aligned approximately perpendicular to the plane of the sky into three basic 
types. The polarization vectors are parallel to filaments when B z ,s/B,p,s ~ 0.1, where B z ,s and 
B<f, t s are respectively the poloidal and toroidal magnetic field components at the outer surface 
of the filament. The polarization vectors are perpendicular to filaments when B z .s / -B<£,s ~ 0.33. 
Intermediate cases result in polarization patterns that contain 90° flips in the orientation of 
the polarization vectors. The flips are symmetric about the central axis for filaments oriented 
parallel to the plane of the sky, but more complicated asymmetric patters result from filaments 
that are inclined at some angle. 



Subject headings: ISM: magnetic fields - polarization - ISM: clouds - MHD 



1. Introduction 

Molecular clouds are highly filamentary structures that are supported against rapid, global collapse 
by a combination of ordered magnetic fields as well as by non-thermal, hydromagnetic turbulence of some 
kind. However, determining the structure of the field has proven difficult. Direct Zeeman measurements 
have been limited to a few points per cloud rather than full-scale maps, and rarely trace the highest density 
gas (Crutcher et al. 1993, 1996, 1999). Polarization maps in optical and near-infrared extinction fail to 
reveal the structure of the field in the dense gas, perhaps as a result of the poor polarizing power of grains 
at these wavelengths in regions of high optical depth (Goodman et al. 1995). 

Observations at far-infrared and sub-millimetre wavelengths have demonstrated that the thermal 
emission from dust grains is often partially polarized (eg. Schleuning 1998, Matthews & Wilson 2000, 
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Hildebrand et al. 1999; see also review by Weintraub, Goodman, & Akeson 2000), which is usually 
attributed to the partial alignment of the emitting grains by the magnetic field. The mechanisms that are 
responsible for grain alignment are not entirely understood, but it is most often argued that dissipative 
processes gradually cause rapidly spinning grains to relax to a state in which grains are, on average, aligned 
perpendicular to the magnetic field (cf. Lazarian, Goodman, & Myers 1997). The possibility now exists, 
therefore, of using completely sampled JCMT polarimetry maps to place strong constraints on the structure 
of the magnetic fields in molecular clouds. 

Current observations reveal the presence of well-ordered magnetic fields within molecular clouds 
and their clumps and cores. However, the results are quite puzzling in that they often show that the 
percent polarization P decreases toward the highest density and most heavily obscured regions (eg. 
Schleuning 1998, Matthews & Wilson 2000; see also Weintraub et al. 2000 and references therein). This 
so-called "polarization hole" effect appears to be an almost universal feature of polarization maps at these 
wavelengths, although its precise origin has been a mystery. Nevertheless, two possible explanations have 
been proposed. Depolarization at high optical depths has often been attributed to the poor alignment 
of grains in dense regions (Goodman et al. 1995). The second possibility is that the effect might be 
due to disordered, sub-resolution scale variations in the orientation of aligned grains in dense regions, 
whose combined polarized emission would tend to mostly cancel (see review by Weintraub et al. 2000, for 
example). Recent high resolution interferometric polarization maps at millimetre wavelengths by Rao et al. 

(1998) have revealed strong, ordered polarization patterns within the polarization hole found by Schleuning. 
The position angle of the polarization vectors changes rapidly on the plane of the sky, which might be due 
to magnetic fields that are twisted on small scales. 

The KL core is embedded in the "integral-shaped" filament of Orion A (See Johnstone & Bally 1999). 
There is some evidence that the polarization hole in the Orion filament is a global feature of the filament, 
and not localized to the KL region. Part of this filament has recently been mapped in 850 /im polarization 
by Matthews and Wilson (2000), who find that the polarization vectors are remarkably ordered and 
generally aligned with the filament. They also find evidence for a decrease in the polarization percentage 
along the symmetry axis of the filament (See Matthews & Wilson 2000, Figure 2). 

This paper demonstrates that a very simple explanation of the polarization hole effect is possible, 
namely that it arises as a consequence of a more general configuration of the magnetic field within 
filamentary clouds than has previously been contemplated by many authors. We show in this work that 
this effect can be rather naturally explained as arising from grain alignment by large-scale, ordered, helical 
magnetic fields threading filamentary molecular clouds. The Matthews and Wilson (2000) data suggest 
that the magnetic field threading the integral-shaped filament is probably very well-ordered. It is therefore 
conceivable that the polarization hole in this region might be due to a large scale, orderly twist in the 
magnetic field, rather than disordered magnetic fields on small scales. In addition, Vallee and Bastien 

(1999) compared the observed 760 \xm polarization patterns of a number of sources with the qualitative 
patterns expected for several field geometries. For the subset of maps which show an ordered magnetic 
field, they argued that helical fields are most consistent with the data. 

We present here, a very simple model of the polarized thermal emission of aligned grains to calculate 
the polarization patterns that would be expected from our helically magnetized filamentary cloud models 
(Fiege & Pudritz 2000a,b; hereafter FP1 and FP2). Polarimetric maps are difficult to interpret because 
very different three-dimensional field geometries can give rise to qualitatively similar polarization patterns 
(see the discussion in Sections 4.1 and ^ for example). This degeneracy is compounded by the high degree of 



uncertainty regarding the composition and distribution of grains, their polarization cross sections, and their 



- 3 - 



alignment efficiencies. However, we show that under the simplest possible set of assumptions, helical held 
models can produce polarization patterns that are not unlike the observed patterns. Our most important 
result is that a large fraction of our models are depolarized along the symmetry axis of the filamentary 
cloud, as a result of the twisted field lines. Thus helical magnetic fields might provide a very simple, and 
purely geometric, explanation for the polarization hole phenomenon. 

2. Filamentary Clouds and Helical Magnetic Fields 

We first summarize the basic properties of our model of magnetized filamentary clouds (see FP1 and 
FP2 for details). Our models assume that filamentary clouds are self-gravitating, magnetized, truncated 
by an external pressure, and in a state of radial quasi-equilibrium between gravity, magnetic stresses, and 
internal pressure gradients. The ISM provides the external pressure for many filaments, but some are 
embedded in molecular clouds (ie. the integral-shaped filament of Orion A). We assume that the magnetic 
field is cylindrically symmetric, with both a poloidal and toroidal component, so that the field is helical in 
general. Note that the toroidal field component must vanish on the filament axis, where the poloidal field 
component is strongest. As one moves out in radius, the winding of the field on each subsequent cylinder 
increases as the toroidal field strength relative to the poloidal field increases. At large radii, the toroidal 
field decreases, but less rapidly than the poloidal field. 

We constructed a three parameter family of models for helically magnetized filamentary clouds, which 
we constrained with an observational data set collected from the literature. We refer the reader to FP1 for 
the precise definition and a full discussion of the model parameters. Briefly, however, our model involves a 
concentration parameter C, which defines the radius of pressure truncation, and two flux to mass ratios 
T z and T^, for the poloidal and toroidal field components respectively, which are assumed to be constant 
within each filament. The models are characterized by a core radius Tq defined as CT(47rGp c ) -1 / 2 , where 
a is the one dimensional velocity dispersion of the gas (assumed to be constant), G is the gravitational 
constant, and p c is the central density along the axis of the filament. The core radius characterizes the 
approximate radial scale at which the density structure changes from being nearly constant inside of ~ r$, 
to an approximately r~ 2 profile in the less dense outer envelope. The fact that we obtain an approximately 
r~ 2 profile is significant because three recent observations (Alves et al. 1998, Lada, Alves, & Lada 1999, 
and Johnstone & Bally 1999) all find density gradients that are consistent with an r~ 2 profile. Our models 
are therefore in good agreement with the existing data on the radial density structure of filamentary clouds. 

3. A Simple Model for Sub-Millimetre Polarization by Aligned Dust Grains 

The model that we adopt for the sub-millimetre polarization arising from our filamentary cloud models 
most closely resembles the analysis by Wardle and Konigl (1990, hereafter WK90) for the polarized emission 
from the Galactic centre disk. The main difference is that our approach allows us to estimate the magnitude 
of the polarization percentage, whereas WK90 worked with normalized quantities. Another difference is 
that we include the possibility of several grain species, with the assumption that the emitting grains are 
a uniformly mixed population with the number density of each individual species proportional to the gas 
density. In addition, we require each grain species to be of uniform temperature and aligned to the same 
extent throughout the cloud. 

Our assumptions regarding the constancy of grain properties and their alignment efficiencies throughout 
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molecular clouds is undoubtedly a simplification of the true state of the grains in some clouds. As an 
example of the debate concerning grain properties, it has often been argued that the grains within dense 
regions of molecular clouds do not efficiently polarize starlight in extinction because they are either nearly 
spherical or poorly aligned at high optical depths (Goodman et al. 1995). There are several reasons to 
think that this might be the case. For instance, grain growth due to coagulation tends to produce slightly 
larger, and possibly more spherical grains in regions of higher density (Vrba et al. 1993), thus decreasing 
their contribution to the polarized emission. This effect might be compounded by the growth of ice 
mantles at high optical depths (Eiroa & Hodapp 1989, Goodman et al. 1995). The nature and efficiency of 
grain alignment mechanisms at high optical depths is perhaps even less certain than the grain properties 
themselves. The alignment might be poor because all efficient alignment mechanisms require supra-thermal 
rotational velocities, and neither the Purcell (1979) mechanism nor the radiative alignment mechanism 
(Draine & Weingartner 1996, 1997) are likely to be efficient at high optical depths. [] Nevertheless, simple, 
quantitative estimates of the alignment efficiency and distribution of grain shapes for various physical 
conditions are presently lacking. It is for this reason that we mainly consider the simplest possible model 
of constant magnetic alignment throughout the cloud. We note however that our models can be readily 
extended to explicitly include radiative alignment. 



Grains result in polarized thermal emission if they are aspherical and partially aligned, most probably 
by a magnetic field. We consider the contribution to the polarized emission from several grain species j, 
whose absorption cross sections parallel and perpendicular to the symmetry axis are given by Gjij and 
C±.j. The thermal emission from these grains is polarized if C\\j and C_lj are not equal. It is convenient 
to define the polarization cross section 



(Lee & Draine 1985, hereafter LD85; WK90). 

To calculate the Stokes parameters /, Q, and U, which completely specify the intensity of the emission 
and its state of linear polarization, we sum the contributions arising from each species of dust within the 
cloud. We assume that species j has temperature Tj, polarization cross section C po ij, and number density 
rij = Cjp, where Cj is a constant and p is the total mass density of the cloud. We also assume that the mean 
alignment of the spin angular momentum with the magnetic field is given by 



A full analysis of the polarized emission applicable at all wavelengths would include the effects of 
absorption and scattering. However, both of these effects can be neglected if we restrict our analysis to 
sub-millimetre wavelengths longer than about IOO/xto since the thermal emission is almost always optically 
thin at these long wavelengths (Hildebrand 1983) and scattering is completely insignificant (Novak et al. 
1989). For example, the peak optical depth of even the dense KL core in Orion is only « 0.25 at 100/im 



1 The radiative mechanism might be important within A v 2 of the surface of the cloud or any embedded sources (Draine 
& Weingartner 1996). 



3.1. Analysis 




Cj_,j — Cy (oblate grains) 

l/2(C|ij — C±j) (prolate grains) 



(1) 




(2) 
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(Novak et al. 1989, Schlcuning 1998), and falls off with increasing wavelength as A - ' 3 , where fi « 1.5 — 2 
(Hildebrand 1983, Andre, Ward- Thompson, & Barsony 2000). Therefore, it is appropriate to treat the 
sub-millimetre polarization as arising purely from emission in most circumstances. 

Following DL85, we define the polarization reduction factor $j for dust grain species j as 

<!\ /,',/• , co- 2 -.. (3) 

where Rj is the Rayleigh polarization reduction factor due to imperfect grain alignment, Fj is the 
polarization reduction due to the turbulent component of the magnetic field, and 7 is the angle between 
the plane of the sky and the local direction of the magnetic field. We assume that Rj and Fj are constants 
for each grain species. The combination of variables C po ij<S>j plays the role of the effective polarization 
cross-section, so that the contributions to the Stokes parameters Q and U from grains of species j are given 
by the following: 

Qj = CpoijRjFjCjB^T^Cjq (4) 
Uj = C po i tj RjFjCjB v (Tj)cjU, (5) 

where B v is the Planck function of the grain at temperature Tj , ds is a distance element along the line of 
sight. The quantities q and u are integrals along the line of sight s, which depend on the density structure 
and the geometry of the magnetic field, but not the properties of the grains. We define the x — y plane as 
being parallel to the plane of the sky, with x pointing West and y pointing North. The angle ip is defined 
as the angle (counter-clockwise) between y and the projection of the magnetic field in the plane of the sky. 
Then, q and u are respectively given by 

p cos 2ijj cos 2 7 ds (6) 
p sin 2^ cos 2 7 ds. (7) 

The angle ip varies as a function of s along the line of sight for 3-dimensional field geometries where the 
local direction of the plane-of-sky field component varies along the line of sight. Note that q and u are 
identical in form to the "relative Stokes parameters" defined by WK90. 

The polarization angle x is determined by solving the equations 

cos2 X =^=L=, (8) 

y<r + u 

U 

sin2 X =^== (9) 

Since equations |§| and |9| are multi- valued, we restrict \ to be on the range [0, w). The special case where 
u = and q remains finite is very important to the analysis that follows. For this case, the polarization 
angle is given by x — when q > and x — 7r /2 when q < 0. 

We proceed to calculate the intensity of the emission, which requires the extinction cross-section from 
equation 3.14 of LD85: 

CextJ — C+.j 1 ~ ( 

where C+ j and ctj are defined as 

°+J - n 



cos 2 7 1 



(10) 
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Note that ctj does not vary with position, since we assume that Rj and Fj are constants for each grain 
species. The contribution of grain species j to the intensity / is given by 

I^C+MT^CjiX-ajZi), (12) 

where £ is the surface density and £2 is a related quantity defined by 

£ = p ds 



S 2 = I Pl^-l ) ds. (13) 



2 3 

Note that £ and £2 depend on the density and magnetic structure of the cloud, but not the grain properties 
The polarization percentage is defined by 



V-^fSL. (») 

where Q, U, and / are obtained by summing equations |], [|, and|l^ over grain species j. It is easy to show 
that equation |lj becomes 

where (a) is a weighted mean of ctj defined as follows: 

= £ y + ^) Cj 

Equations ^, 0, and |l^ show that the polarization pattern is determined by the density and magnetic 
structure of the gas, plus a single parameter (a) related to the grain cross-sections and alignment properties. 
We estimate (a) as follows. We assume that the optimal magnetic field geometry, in which 7 = and 
ip = const, gives rise to the maximum polarization percentage p m ax that is normally observed. The 
maximum polarization percentage is obtained from equation ^|, with the help of equations ^ and ^: 

- ^ 

1 - (a)/6 

Since the maximum polarization percentage observed in sub-millimetre polarization maps is rarely greater 
than 10%, equation |l7| suggests that 0.1 is a reasonable choice for (a). Using this choice, our approach 
predicts the magnitude of the polarization percentage obtained from our models, which is generally less 
than Pmax- This is in contrast to WK90 who work with normalized quantities. We caution the reader that 
(a) could, in principle, vary from region to region, which would affect the magnitude of the polarization 
percentage that we predict. It would also slightly affect the shapes of the polarization profiles, since (a) 
apppears in the denominator of equation hq and cannot be factored out. 



4. Results 
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4.1. Transverse Fields Across Filamentary Clouds? 

Polarimetry only maps the component of the magnetic field parallel to the plane of the sky. Suppose 
that the magnetic field threading the integral-shaped filament in Orion A contains only this plane-of-sky 
component. If this were true, the Matthews and Wilson (2000) map, as well as the smaller scale map by 
Schleuning (1998), would suggest that we are seeing either a filament or an edge-on sheet being impaled 
by a well-ordered field that is perpendicular to the filament axis or the midplane of the sheet. While this 
might explain the overall sense of the polarization vectors, it cannot easily account for the depolarization of 
these maps seen toward the centre of the putative sheet or filament. For such a transverse field model, it 
is easy to verify from equations [| ^, and ^| that p would be constant over the entire filament, since C = 



everywhere and ip — const. The analysis of Section 3.1 would not predict the observed depolarization along 



the axis of the filament or the midplane of the sheet. It could, however, be accomplished if the polarization 
were due to grains that are preferentially more spherical or poorly aligned in dense regions. 



We show in Section 4.2 that helical fields threading a filamentary cloud can result in depolarization 
toward the filament axis, otherwise similar in appearance to what would be expected from the above 
transverse field scenario. However, the depolarization is due to the 3-dimensional structure of the field in 
this case and does not depend on the grain shapes or their alignment. 



4.2. Polarization Patterns and Depolarization Effects for Helical Fields 

We present polarization maps of our helically magnetized models of filamentary clouds using the 



method discussed in Section (SdJ. We first restrict our parameter space by only showing results for filaments 



whose axes lie in the plane of the sky. In Section 4.4, we show maps for a filament at several inclination 
angles on the sky. 

Figure [l] shows three representative models to illustrate the general types of behaviour that we find 
in our maps. We refer to these qualitative patterns as types 1 to 3 for simplicity, although we emphasize 
that the underlying models represent a continuum in parameter space (see FP1). In panel a) we show a 
model (type 1) where the polarization vectors are everywhere parallel to the filament. The most striking 
feature of the map is the depolarization along the axis. Qualitatively, this map shares some features with 
the Matthews and Wilson (2000) map of the Orion filament, specifically the overall orientation of the 
polarization vectors and the depolarization toward the central regions. Panel b) shows a model (type 
2) whose polarization vectors are oriented opposite to the type 1 model in panel a). If such a pattern 
were observed, it could be misinterpreted as the result of a purely poloidal field. These models have 
two depolarized regions, with one on either side of the filament axis. The polarization on the filament 
axis is a local maximum, although the strongest polarization in the map is usually at the surface of the 
filament. Panel c) shows a more complicated case (type 3), where there are 90° flips in the orientation of 
the polarization vectors. The polarization percentage passes through zero at all locations where a flip in 
orientation occurs. We classify any patterns containing such flips in orientation as type 3. 

It may seem surprising that filamentary clouds that are threaded by helical fields always have 
polarization vectors that are aligned parallel or perpendicular to the filament. This result is obtained 
regardless of whether the filament is parallel to the plane of the sky or inclined at at some angle. It is 
easy to understand this by the following argument (see Carlqvist & Kristen (1997) and FP1). Any line of 
sight passing through a filament encounters a given magnetic flux tube twice; first on the front side of the 
filament, and then on the back. By symmetry, the plane-of-sky components of the magnetic field vectors 
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are mirror images of each other at these positions, so that if) is an odd function of the position s along 
each line of sight. Therefore, equations |^ and [7] show that the contributions to u exactly cancel, while the 
contributions to q add in equal proportion. 

From the definition of the polarization angle (equations || and ^), we observe that the polarization 
vectors are oriented perpendicular to the filament (x = 0) when q > and parallel to the filament (\ — 7r/2) 
when q < 0. The dominant magnetic field component on each line of sight determines the orientation of the 
field vectors. This is easily seen by the following argument. Consider first a purely poloidal field threading 
a filament oriented in the North-South direction, parallel to the plane of the sky. The Stokes vector q > 
in this case, so that the polarization vectors are perpendicular to the filament, since cos 2?/; = 1 everywhere. 
On the other hand, if the field were purely toroidal, cos 2-0 = —1, so that q is negative and the polarization 
vectors are parallel to the filament. Now consider the case of a helical field threading a filament oriented 
parallel to the plane of the sky. It is obvious that the poloidal and toroidal components of the magnetic 
field provide competing contributions to q through 0, and that the polarization direction at each position 
is determined by the field component that makes the dominant contribution to equation || along each line 
of sight. 

The competing contributions to q from the poloidal and toroidal field components may partially or 
completely cancel at some positions, whenever the poloidal and toroidal fields are of comparable strengths. 
This is the explanation for the depolarized regions seen in all three panels of Figure [j]. 



4.3. Exploration of the Parameter Space 

We explore our parameter space in the same spirit as we did in FP1, FP2, and Fiege & Pudritz 
2000c (Hereafter FP3). We compute polarization maps for models scattered randomly throughout our 
3-dimensional parameter space (C, T z , T^), subject to the observational constraints discussed in FP1. 
A slight complication is that we exclude some models that were found to be highly unstable to short 
wavelength, rapidly growing MHD instabilities in FP2. We note that a similar criterion was used in 
FP3 to select "parent" filaments for our prolate core models. Specifically, we exclude models for which 
~ UJ maxl i^Gpc) > 0.05, where — ojf nax is the maximum squared growth rate, and p c is the central density 
of the filament (See FP2, Figure 13). 

In Figure||, we show how the polarization pattern depends on B z s / B(f>,S, where B z ^s and -B^s are 
respectively the poloidal and toroidal field components, evaluated at the surface of the filament: 



B z ,s _ e-°T 



B. 



<j>,s 



(18) 



The dots, squares, and x's represent models of types 1, 2, and 3 respectively. Note that -B^s/i^s is small 
for many of our models because the toroidal magnetic field dominates in the outer envelope, near the radius 
of pressure truncation. The ratio of B z /B^ is substantially higher in the interior regions of the filament. 

An interesting feature of Figure || is that the polarization patterns are of the first type for a large 
portion of our models. Thus, many of our models are qualitatively similar to the Matthews and Wilson 



(2000) map of the Orion filament, discussed in Section 4.2. We find this type of pattern for most models 



with B z ^s/B^,s ~ 0.1. Generally, polarization patterns of the second type occur when B Zi s/B^^s ~ 0.33, 
and the third type occurs for intermediate values between about 0.1 and 0.33. This may be understood as 
follows. When Bzs/B^s is small, the models are dominated by the toroidal field component so that the 
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poloidal field is ineffective at canceling the polarization due to the toroidal field, as discussed in Section 



4.2. Thus, the polarization pattern is of type 1. However, as B z g is increased relative to B^.g, the 
competition becomes stronger until B z g becomes dominant along some line of sight, which first occurs 
when B z g/B^g ps 0.1. This results in a 90° flip in the polarization vectors at this position, which we 
categorize as the first Type 3 model. The exact opposite occurs when B z .g / B^^g is increased past about 
0.33. Past this point, the toroidal field component becomes too weak compared to the poloidal field to 
produce a flip in the polarization vectors resulting in Type 2 patterns. 

Figure [| shows the degree to which the emission is depolarized for all three types of polarization 
pattern. We define the polarization hole depth as (P ma x — Pmin) /Pmax, where P m ax is the maximum 
polarization percentage in the map and Pmin is the local minimum polarization percentage at the location 
of the polarization hole with the lowest polarization. We find that the polarization hole depth generally 
increases, with scatter, as a function of B Zt g/Bj, g for Type 1 models from 0% when B z g/B^ g = to 100% 
when Bzs/B^g = 0.1. This is easily understood by essentially the same argument given in the previous 
paragraph. There is no significant polarization hole when B z .g / B^^g is small because contributions to 
the polarization arising from the poloidal field do not effectively cancel the larger contributions from the 
toroidal field along any line of sight. The depth of the polarization hole increases until B z ^gjB^,g w 0.1, 
where the first Type 3 pattern emerges. Note that the polarization hole depth is always 100% for Type 
3 models, since P m in = at the locations where the oriention of the polarization vectors flips by 90°. 
Increasing Bz^g/B^^g decreases the depth of the polarization hole for Type 2 models, since the contribution 
to the polarization from the toroidal field becomes progressively less effective at canceling the polarization 
due to the dominant poloidal field. 

We define the width of the depolarized region as the distance between the polarization maxima on 
either side of the polarization hole. Panel [|b shows the ratio of this width divided by the filament diameter. 
Generally, this ratio increases with B Zt g/B^ } g for Type 1 models, from nearly to about 0.4, which occurs 
at the point where the models change to Type 3. This fractional width continues to increase past this 
point until B Zt g / B^^g k, 0.23, where the polarization width jumps discontinuously to the full width of the 
filament. This happens because of a change in the qualitative behaviour of the polarization vectors near 
the edge of the filament. Note that the polarization percentage shown in Figure ^ is a minimum at the 
outer edge of the filament for Type 1 models, but a maximum for Type 2 models. As we move through a 
sequence of Type 3 models, the polarization at the edge gradually increases as a function of B z ^gjB$^g until 
B Zl g/B c f >t g ks 0.23. Past this point, the maximum polarization for the filament occurs at the edge, so the 
width of the depolarized region becomes equal to the width of the filament. 



4.4. The Effect of Inclination Angle 

The polarization patterns due to helically magnetized filaments become more complicated for filaments 
that are inclined relative to the plane of the sky. As an example, in Figure |4| we show maps of one filament 
tilted 10°, 30°, and 60° relative to the plane of the sky. The filament model shown in Figure ^ is the 
same one as shown in the top two panels of Figure [l]. Each polarization vector always remains parallel or 
perpendicular to the symmetry axis for the reasons discussed in Section [T| above, but the patterns become 
asymmetric relative to the central axis of the filament. It is particularly notable that filaments inclined 
relative to the line of sight usually result in asymmetric polarization patterns that include a 90° flip in the 
orientation of the polarization vectors on one or both sides of the filament. 
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5. Discussion Sz Summary 

The main purpose of this paper is to calculate polarization maps for the models of filamentary 
molecular clouds that we presented in FP1 and FP2. When the filament is oriented parallel to the plane of 
the sky, the polarization patterns resulting from our models can be qualitatively classified into three main 
types of behaviour. The polarization vectors are parallel to the symmetry axis of the filament (Type 1) 
when the toroidal field is dominant, and perpendicular (Type 2) for models whose poloidal field dominates. 
The third type of pattern is more complicated, with polarization vectors that flip from being aligned parallel 
to the filament axis to perpendicular at some radius or radii (in projection). Generally, the first type of 
pattern occurs when B Zi s/B^, t s ~ 0.1, while the second type occurs 

B z ,s/B^ s ~ 0.33. The third type of 
polarization pattern occurs at intermediate values of B Zj s/B^, t s Our most important result is that helical 
fields result in "polarization holes," in which the emission is depolarized at some positions in the interior of 
the filament. Our Type 1 models qualitatively agree with the polarization structure of the "integral-shaped" 
filament in Orion A (Matthews & Wilson 2000). 

Many of our models result in polarization patterns that contain 90° flips in the orientation of the 
polarization vectors. The flips are symmetric about the symmetry axis when a filament with comparable 
poloidal and toroidal magnetic field strengths is oriented parallel to the plane of the sky. Filaments that are 
inclined relative to the line of sight generally result in asymmetric polarization flips on one or both sides of 
the central axis of the filament. 

We have considered only models of non- fragmented, cylindrically symmetric filaments in this paper. 
This is obviously an idealization of real filamentary clouds, many of which have suffered gravitational 
fragmentation and formed strings of embedded cores (cf. Schneider & Elmegreen 1979, Dutrcy ct al. 1991; 
see also FP2 and FP3), and may also deviate from perfect cylindrical symmetry. We will address the 
polarization patterns of embedded cores in a future analysis. It is difficult to comment on the general effects 
of asymmetry. However, it is worthwhile to consider the simplest type of non-axisymmetric perturbation, 
that corresponding to a kink mode in which the filament is bent into a transverse sinusoidal wave. If the 
sinusoid lies in the plane of the sky, the front to back symmetry is preserved and the polarization vectors 
remain locally parallel or perpendicular to the filament. The depolarization along the central axis of the 
filament relies on a cancellation between contributions to the polarization from the central "backbone" 
of poloidal flux and the surrounding envelope dominated by the toroidal field component. Since these 
features would be preserved in such a perturbation, the polarization hole would also likely remain. If, on 
the other hand, the sinusoidal perturbation is perpendicular to the plane of the sky, then the front to back 
symmetry is broken and the polarization vectors would no longer be strictly parallel or perpendicular to 
the filament. A segment which is bent away from the observer is more compressed on the front side and 
more rarified on the back, so the front side would contribute more to the polarization. Assuming that the 
filament is oriented in the north-south direction and the field lines have right-handed (positive) helicity, the 
polarization vectors would shift toward the north-east/south-west plane if they were originally oriented in 
the north-south direction. The polarization hole would be preserved for the same reasons as in the case of a 
perturbation parallel to the plane of the sky. Thus, we see that non-axisymmetric distortions of a filament 
containing a helical field still ought to show a "polarization hole" . 

Ignoring systematic variations of grain shapes and alignment efficiencies with density, we have shown 
that the "polarization hole" effect has a very simple and entirely geometric explanation if the field lines are 
twisted into helices. It might alternately be explained if the emitting dust grains are poorly aligned in the 
central regions of filamentary clouds, but this requires a greater understanding of dust grain physics. 
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More rigorous tests of our models will require a full set of polarization maps of filamentary clouds. 
Such a data set, when combined with well determined radial density profiles, will provide an excellent 
challenge to all cloud models, not only our own. Readers wishing to use our models to compare with data 
should contact J. Fiege to obtain fits of our model to data, or polarization maps for any inclination angle 
and set of input parameters (C, T z , and T^). It is probably true that even with the density profiles and 
polarization structure of filamentary clouds determined by observations, it would still be necessary to make 
spatially resolved Zeeman measurements in order to conclusively prove the existence of helical fields. By 
making such maps, we might one day decide whether grain properties or the structure of the magnetic field 
is responsible for the polarization hole effect. 
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x/r Q x/r Q 

Fig. 1. — We show examples of the three types of polarization maps that we obtain for filaments oriented 
perpendicular to the line of sight. (Top row:) Type 1: A toroidally dominated field pattern (C = 0.62, 
T z = 5.1, T = 13.6). (Centre row:) Type 2: A poloidally dominated field pattern (C = 0.31, T z = 8.7, 
r = 11.9). (Bottom row:) Type 3: A mixed field pattern (C = 0.54, T z = 10.4, T = 11.4). 
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Fig. 2. — We show the polarization percentage as a function of B Zi s/B^s- According to the classification 
scheme discussed in Section 4.2, dots represent type 1 patterns, squares represent type 2 patterns, and x's 
represent type 3 patterns. 
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Fig. 3. — We show the relative depolarization of our models (left panel), and the width of the depolarized 
region (right panel). The dots, squares, and x's use the same classification scheme as in Figure |[ 




Fig. 4. — We show the polarization patterns expected for one of our models oriented at 10° (top row), 30° 
(centre row) , and 60° (bottom row) relative to the plane of the sky. The model used here is the same one 
used in the top row of Figure pi 



